Cloning of Gaussian states by linear optics 
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We analyze in details a scheme for cloning of Gaussian states based on linear optical components 
and homodyne detection recently demonstrated by U. L. Andersen et al. [Phys. Rev. Lett. 94, 240503 
(2005)]. The input-output fidelity is evaluated for a generic (pure or mixed) Gaussian state taking 
into account the effect of non-unit quantum efficiency and unbalanced mode-mixing. In addition, 
since in most quantum information protocols the covariance matrix of the set of input states is not 
perfectly known, we evaluate the average cloning fidelity for classes of Gaussian states with the 
degree of squeezing and the number of thermal photons being only partially known. 
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I. INTRODUCTION 



The generation of perfect copies of an unknown quan- 
tum state is impossible according to the very nature of 
quantum mechanics. This is succinctly formulated by the 
no-cloning theorem It is, however, possible to 

make approximate copies of a quantum state by using a 
quantum cloning machine . Originally, such a machine 
was proposed for cloning of qubits and has later been 
demonstrated experimentally 'u\. Shortly after this de- 
velopment, a continuous variable (CV) 7J analog of the 
qubit quantum doner was proposed [3, Q and recently 
it was shown that a CV optimal Gaussian doner of co- 
herent states can be implemented using an appropriate 
combination of beam splitters and a single phase insensi- 
tive parametric amplifier [lol[T]| . Although this proposal 
sounds experimentally promising, the implementation of 
an efficient phase insensitive amplifier operating at the 
fundamental limit is a challenging task. This problem 
was solved by Andersen et al. who proposed and 
experimentally realized a much simpler configuration for 
optimal cloning of coherent states. The realization relies 
on simple linear optical components and a feed-forward 
loop. As a consequence of the simplicity, as well as the 
high quality of the optical devices used in this experi- 
ment, performances close to optimal ones were attained. 
In turn, the resulting cloning machine represents a highly 
versatile tool for further investigations on transformation 
of quantum information from a single system to many 
systems. 

A commonly used figure of merit to quantify the per- 
formance of cloning machines is the fidelity which is a 
measure of similarity between the hypothetically perfect 
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clone, i.e the input state, and the actual clone. If the 
cloning fidelity is independent on the initial state the ma- 
chine is referred to as a universal doner. On the other 
hand, if the efficiency of the cloning action depends on the 
input state, then the proper measure in order to assess 
the performances of the machine is the average fidelity, 
which weight the fidelities associated to possible input 
states with the corresponding occurrence probability. In 
other words, for non-universal doners, the alphabet of 
input states, and the distribution thereof, must be taken 
into account while evaluating the fidelity Such an aver- 
age fidelity has been considered in 0,0; E|- However, 
in all these references it is assumed that the input alpha- 
bet is only consisting of coherent states, hereby keeping 
the covariance matrix of all the possible input states con- 
stant. On the other hand, in some experimental realiza- 
tions, the covariance matrix is not perfectly known due to 
uncontrollable fluctuations, and therefore it is important 
to include this uncertainty into the analysis. 

The aim of this paper is two-fold. At first we present a 
thorough theoretical description of the cloning machine 
described in Ref. using a suitable phase-space anal- 
ysis. In this way the full quantum dynamics of the ma- 
chine can be taken into account; in particular we include 
the effect of losses in the detection scheme, as well as 
variations in the setups beam splitter ratios. The second 
topic of the paper is to investigate the average fidelity 
of the cloning machine for different ensembles of input 
states such as sets made of displaced squeezed or dis- 
placed thermal states with the squeezing parameter, or 
the number of thermal photons, distributed according to 
predefined distributions. 

The paper is structured as follows: in Sec. we re- 
view the main components of the cloning machine based 
on linear optics, whereas in Sec. IIIII we calculate the 
input-output fidelities for the case of generic Gaussian 
states, and for specific classes including coherent, dis- 
placed squeezed and displaced thermal states. Finally, 
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Sec. IIVI closes the paper with some concluding remarks. 



II. THE LINEAR CLONING MACHINE 



Optimal Gaussian cloning can be realized us ing a phase 
insensitive amplifier and a beam splitter pil ml. How- 
ever, it has been recently shown, theoretically and exper- 
imentally, that the parametric amplifier can be replaced 
by a simpler scheme involving only linear optical compo- 
nents, homodyne detection and a feed- forward loop |l2l |. 
This scheme, which is schematically depicted in Fig. ^ 
will be referred to as the linear cloning machine through- 
out the paper. 




FIG. 1; Cloning of Gaussian states by linear optics: the in- 
put state Qin is mixed with the vacuum go at a beam splitter 
(BS) of transmissivity ri. The reflected beam is measured 
by double-homodyne detection and the outcome of the mea- 
surement a; -I- zy is forwarded to a modulator, which imposes 
a displacement g{x + iy) on the transmitted beam, g being 
a suitable amplification factor. Finally, the displaced state 
is impinged onto a second beam splitter of transmissivity T2. 
The two outputs, qi and Q2, from the beam splitter represents 
the two clones. 

The input state, denoted by the density operator 
Pin, is mixed with the vacuum at a beam splitter (BS) 
with transmittivity ti. On the reflected part, double- 
homodyne detection is performed using two detectors 
with equal quantum efficiencies 77: this measurement is 
executed by splitting the state at a balanced beam split- 
ter and, then, measuring the two conjugate quadratures 
X = -^{a + a)) and y = "^(Q'"^^), with a and being 
the field annihilation and creation operator. The out- 
come of the double-homodyne detector gives the complex 
number a — x + iy. According to these outcomes, the 
transmitted part of the input state undergoes a displace- 
ment by an amount ga, where g is a suitable electronic 
amplification factor, and, finally, the two output states, 
denoted by the density operators qi and Q2 , are obtained 
by dividing the displaced state using another beam split- 
ter with transmittivity T2. When ti — T2 — 1/2 , g — 1 
and ?7 = 1, the scheme reduces to that of Ref. [ij, which 
was shown to be optimal for Gaussian cloning of coher- 
ent states on the basis of a description in the Heisenberg 
picture. Here we apply a different approach which cap- 
tures all the essential features of the machine. Towards 



this aim, in the following we carry out a thorough de- 
scription of the machine using the characteristic function 
approach. 

The characteristic function Xm 

(Ai) = xiftnKAi) asso- 
ciated with a generic Gaussian input state gin of mode 1 
reads: 



Xi: 



,(Ai) = exp|-iAfcrin Ai - iAfXinj , (1) 



where Ai — (xi,yi)-^, (• • • )^ denotes the transposition 
operation, and 



711 712 
721 722 



(2) 



with 712 — 721, is the covariance matrix. Xin = 
Tr[gi„ is the vector of mean values, x and y be- 

ing the quadrature operators defined above. The vacuum 
state po — |0) (0| of mode 2 is described by the (Gaussian) 
characteristic function 



Xo 



(A2) =xbo](A2) = cxp{-iA^croA2} , (3) 



where ctq ~ 5I2, I2 being the 2x2 identity matrix. In 
turn, the initial two-mode state Q = Qm® Qo is Gaussian 
and its two-mode characteristic function reads: 



x[^)](A) = exp{-iA^^A-^A^x} , 



with 



Cin 








To 



X 



(Xi„,0) 



(4) 



(5) 



and A = (Ai, A2)"^. Under the action of the first BS the 
state x[p](A) preserves its Gaussian form, namely 

x[e](A)^xb'](A) = exp{-iA^^A-7A^x} , (6) 

where g' = Ubs,i Qin 'S) Qo U^g ^, while its covariance ma- 
trix and mean values transform as [T^ : 



A 


c 




B 



X X = Sgg 1 X = (Xl, X2) , 

A, B, and C are 2x2 matrices, and 



'S'bs, 



Vnl2 


Vl - n I2 


-Vl - n I2 


^/^l2 



(7) 
(8) 

(9) 



is the symplectic transformation associated with the evo- 
lution operator J7bs,i of the BS with transmission t\. 
Note that (J is an entangled state if the set of states to 
be cloned consists of nonclassical states, i.e. states with 
sing ular Glauber P-function or negative Wigner function 

III El. 
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The subsequent step is to describe double-homodyne 
detection with quantum efficiency rj on the reflected 
beam. This action can be described by the following 
positive operator- valued measure (POVM): 

where = {l — r])/ri and |^) is a coherent state. Eq. itTUI) 
describes a Gaussian measurement, the characteristic 
fmiction associated with 11,, (a) has the form 



^[n,,(a)](A2) = -exp|-iA^crM A2 - iA^ Xm^ , 
with Xm — (Re[a], Im[a]) and 



1 



2-77 



(12) 



2 " 2?7 

The probability of obtaining the outcome a is given by 
-Tri2[f?'l(»n,,(a)] (13) 
d4Ax[f>'](A)x[II®n,(a)](-A) (14) 



1 



(27r)^ 7r4 

_ exp {-^(Xm - X2V (Xm - X^)} 
7r^Det[S] 

(15) 

where x[I ® n,,(a)](A) = x[I](Ai) x[n,,(a)](A2), 
x[I](Ai) = 27r(5(2)(A^) and (^^^^C) is the complex Dirac's 
delta function. We also introduced the 2x2 matrix 
S = B + (Tm- 

The conditional state Qc of the transmitted beam, ob- 
tained when the outcome of the measurement is a, i.e., 



Tr2[g'n„(a)] 



(16) 



has the following characteristic function (for the sake of 
clarity we explicitly write the dependence on Ai and A2) 



x[ec](Ai) 



^2^^^ x[g1(Ai,A2)x[n„(a)](-A2) ^^^^ 



: exp 



|-iAf A-CS-^C^ Ai 
iX^ X2 + lA^ [CS-i X2 - Xi] } 
exp{-iX^S-iXM 

}. (18) 



iS-iX2 



S^^C^Ai 



Now, the conditional state Qc is displaced by the amount 
ga resulting from the measurement amplified by a fac- 
tor g. By averaging over all possible outcomes of the 
double-homodyne detection, we obtain the following out- 
put state: 



£»d = / (fapr,{a)D{ga)QcD\ga) 
Jc 



(19) 



with D(z) being the displacement operator. In turn, the 
characteristic function reads as follows: 

x[ed](Ai) = 2exp{-iAf<TdAi-zAfXd} , (20) 

with (Td^ A + gCE + 2C^) and Xd= Xi+ 5X2. The 
conditioned state p9|) is then sent to a second beam split- 
ter with transmission T2 (see Fig. where it is mixed 
with the vacuum qq, and finally the two clones are gen- 
erated. Note that, in practice, the average over all the 
possible outcomes a in Eq. ((T^ should be performed at 
this stage, that is after the second beam splitter. On the 
other hand, because of the linearity of the integration, 
the results are identical, but performing the averaging 
just before the beam splitter simplifies the calculations. 
Since is still Gaussian, the two- mode state g{ = Qd'^go 
is a Gaussian with covariance matrix and mean given by 



(Tf = 



CTd 








era 



Xf = (Xd,0)^, (21) 



respectively, which, as in the case of Eqs. Q and lO, 
under the action of the BS transform as follows: 



(Tf (Tout = 5^33,2 T SbS,2 



Xf Xout — S^BS 2 — ■^2)'^ , 









^2 



(22) 
(23) 



where £/k and are 2x2 matrices, and 5^68.2 is the 
symplectic matrix given by Eq. @ with ti replaced by 
T2. Finally, the (Gaussian) characteristic function of the 
clone pfc, /c = 1, 2, is obtained by integrating over A^, h ^ 
k, the two-mode characteristic function x[Pout](Ai, A2), 

where gout = Ubs.2 Qf ® Qo C^bs,2' i-^-: 

xH(Afe) = / d^Ah xkut](Ai, A2) (24) 
= exp{-iAf j24Afe -zA^^fc} . (25) 

Let us now focus our attention on Xout^ the explicit 
expressions of and are 

^1 = V^(\/^ + 5\/T^)^in, (26) 

J^2 = \/r^(Vn+.g\/T^n)Xi„. (27) 

As a matter of fact, in order to have two output Gaussian 
states with the same means = one should put 
T2 — 1/2; furthermore, if one also sets 
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1-Tl 



l-Tl 



(28) 



then = 3^2 = X[^, corresponding to unity gain 
cloning. On the other hand, £/k can be written in a 
compact form as follows: 

M = 5(1-^2) l + r2r((Tin), (29a) 

^2 = ir2 1 + (1 - T2) r((Tin) , (29b) 
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where 



with 



^(7) = l-ri+g Ti-2\/(l-ri)Ti+A2 +^(7) 



^(7) = [n + .g (1 - ri + 2v/(l-ri) 



(31) 
(32) 



Now, if T2 — 1/2 and g — gs, one has = ^2 and ,^1 = 
as we have seen above, i.e., the cloning becomes 
symmetric. Furthermore, when also ri = 1/2, thanks to 
Eqs. (|25|l and (|29|l we have that the cloning map for the 
scheme in Fig.^is given by the following Gaussian map: 



,(£'in) 



exp 



^1 D{^)g,,,D\^)., 



'GN 



(33) 



where = _| + . Finally, although does not 

appear in Eq. H25|) . for the sake of completeness, we give 



its analytic expression: 



- v/(l - ^2)t2 [il-r(cri, 



(34) 



In the following we will analyze the input-output fideli- 
ties for a generic (pure or mixed) Gaussian state. In par- 
ticular, we will consider three classes of Gaussian states, 
i.e. coherent, displaced squeezed and displaced thermal 
states. 



III. CLONING OF GAUSSIAN STATES 
A. Fidelity 

Usually, the performance of cloning machines are quan- 
tified by the fidelity which is a measure of the similarity 
between the hypothetically perfect clone and the actual 
clone. In its most general form, the fidelity is given by 
the Uhlmann's transition probability 0| 



F{Qin, Qk) = Tr 



'Qin Qk \/ Qin 



(35) 



and satisfies the natural axioms 

• F[Qin,Qk) < 1 and F{gin,gk) = 1 if and only if 

Qin — Qk', 

• F{gin, Qk) = F{gk, ftn); 

• if Qin is a pure state Qm — {ipinjiipinl, then we have 

F{gin, Qk) = (V'inlfi'fclV'in); 

• F{gin, gk) is invariant under unitary transforma- 
tions on the state space. 



Furthermore, when g^ and g^ are Gaussian states of the 
form (HJ and the fidelity ^ becomes jMEJl 

F^ = F{gin, gk) = , , , 7? 

X exp {-i(Xin - ^fc)^(cri„ -f ^fe)-i(Xin -^k)] , 

(36) 

where 5 = 4(Det[<Tiii] - i)(Det[.e4] - i). Note that for 
pure Gaussian states Det[(Tin] = i, and in turn 6 = 0. 



For symmetric cloning, i.e. for T2 
Eq. then Eq. ^ reduces to 



F,= 



1/2 and 3 = <?s in 



(37) 



In general, the cloning fidelity in 1)36(1 is state depen- 
dent, and therefore the figure of merit to be considered 
is the mean cloning fidelity, averaged over the ensemble 
of possible input states. In order to evaluate this quan- 
tity, we parametrize the input ensemble (class) {^?in(A)} 
of different Gaussian states, by A G and consider each 
of them occurring with the a priori probability p( A) . The 
average fidelity then reads 



f dXp{X)F^{X). 
Jn 



(38) 



Within the set of possible states both the mean values 
as well as the covariance matrices may vary. Assuming 
that the probability distribution p{\) is factorisable into 
a distribution for the mean values p{a) and a distribution 
for the covariance matrix, p(<Tin), we may write p(A) = 
p{a) p{crin), and the average fidelity reads 



Fr, = d(Tindap{a)p{(Tin) Fjj{crin,a) . 

Jn 



(39) 



In the extreme case where both (Tin and a are fixed, the 
input state is completely known and perfect cloning with 
unit fidelity is of course possible. A more interesting 
scenario is when the covariance matrix is fixed, as for 
example the case in which the set is made by coherent 
states, while the displacement (that is, the mean value) 
is random. In this case the average fidelity reduces to 



Fn^ 



d'^ap{a) Fjj{a) 



(40) 



If Ti = 1/2 and g = gs, the map (|33|l is covariant 
with respect to displacements, meaning that if two in- 
put states are identical up to a displacement their re- 
spective clones should be identical up to the same dis- 
placement ,22iJ. Indeed, if the input state is of the form 
Qinict) = D{a) gsD^a), gs being a seed state, then the 
fidelity Fr^{a) actually does not depend on complex pa- 
rameter a, and, as consequence, Fn = Frj. Therefore, 
in this case the noise added by the cloning process 
is the same noise added in cloning of coherent states, 
i.e., the cloning is optimal. Notice that the correspond- 
ing optimal fidelity F is not necessarily equal to 2/3 [see 
Eq. 
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B. Coherent states 

Before addressing the general case let us reconsider 
cloning of (pure) coherent states. For this set of states 
our linear machine provides universal cloning, i.e. state 
independent fidelity. In Fig.|5]we plot the fidelity as given 
by Eq. (|37|l . as a function of ti, for different values of rj 
and for T2 = 1/2, g = gs. In this case, corresponding 
to symmetric cloning, the machine yields the optimal fi- 
delity F = 2/3 predicted for universal Gaussian cloning 
of coherent states. Notice that the optimal fidelity is 
achieved with ri = 1/2 and rj — 1; by expanding the fi- 
delity up to the second order around ri = 1/2 we obtain 



F = F 



2ri 



1-K277 



1 - 



1 + 2t] 



^From this expression we clearly see that the cloning ma- 
chine proposed in is robust against fluctuations of the 
BS ratio. This conclusion can be also directly drawn from 
Fig.El 



We have now seen that by fixing the covariance ma- 
trix of the input states to coherent states, the fidelity is 
a function of the distribution (being delta, flat or Gaus- 
sian) of these states. This aspect has been investigated 
in the literature |l3|- In contrast, the case where the co- 
variance matrix may fluctuate has not received much at- 
tention heretofore. In the following sections we therefore 
discuss the average cloning fidelity for classes of states 
with covariance matrix randomly distributed according 
to a predetermined distribution. We assume that the 
displacement of the input state is random and that the 
doner is set to unity gain (that is invariant with respect 
to the displacement corresponding to g = gs)- In this 
case the average over the mean value is trivial and the 
average fidelity can be written as 



Fr, = I dcTu 



(44) 



C. Squeezed states 



1 
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FIG. 2: Linear cloning fidelity Fr, for a coherent state as a 
function of the BS transmittivity n for different values of the 
quantum efficiency rj: from top to bottom rj = 1.0, 0.75, and 
0.5. We set T2 = 1/2 and g ~ Qa (symmetric cloning). The 
dashed line is the value 2/3. The fidelity does not depend on 
the coherent state amplitude. 

Let us however note that the fidelity value F = 2/3 is 
the optimal one only if the input distribution of coherent 
states is flat, that is, if we have no a priori information 
about the amplitudes. If the set of coherent states is 
restricted such that the distribution of amplitudes is the 
Gaussian 

Pa(«) = ^exp(-^l , (41) 

the average fidelity can be increased by choosing a dif- 
ferent gain However, in this scenario the cloning 
action becomes state dependent, and the integration in 
(|40|l should be explicitly performed. By optimizing the 
gain we find ^3 

^=TT^' -^>l + V2, (42) 

F= — , if al<l + V2. (43) 

2 + (3-2V2)a2' " ^ ' 



When the input Gaussian state is the squeezed state 
\a,£,) = D{a)S{£,)\0), where £> (a) = exp{aa''' - a*a} and 

S'(^) — exp{i(^a'f^ — ^*a^)} are the displacement and 
squeezing operator, respectively, the entries of the input 
covariance matrix (j^J are 

711 = i (cosh2|^| -f sinh2|^| cos(p) , (45a) 
722 = 5 (cosh2|^| - sinh2|e| cos(^) , (45b) 

712 = 721 = -5sinh2|^| sin.^, (45c) 

where we put ^ = l^le**^; obviously, when ^ = the 
squeezed state |a,^) reduces to the coherent state \a) 
and <Ti,i = ^l2- Note that in this Section we are ad- 
dressing the case of an unknown squeezing parameter ^ 
(randomly distributed according to a given probability 
density): when it is known, the optimal strategy in the 
Gaussian regime is to perform the unsqueezing operation 
5^^(^) just before the cloning machine, proceed as in the 
case of coherent states, and, at the output stage, apply 
the squeezing operation S{^) to both the clones which 
yields a fidelity of 2/3 (independent on the amount of 
fixed squeezing) as in the coherent state case 10]. 

However in the case of an unknown squeezing param- 
eter the squeezing action 5'(^) is not known. Therefore 
in the following, we investigate the cloning of unknown 
squeezed states using the cloning machine outlined in this 
paper. First we note that since the linear elements in- 
volved in the cloning machine do not affect the phase 
of the input state, the fidelity i^^(^) depends only on 
1^1 and, without loss of generality, we may take ^ as real. 
The fidelity for Gaussian squeezed input states, using the 
coherent state cloning machine, is given by 

m) = , ^ ■ (46) 

J(5 + 2A2)2 -I- 16(1 + 2A2) sinh^ |^| 
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FIG. 3: Plot of the fidelity of the squeezed state \a, ^) 

as a function of |^| for different values of rj: from top to bottom 
7] = 1.0, 0.75, and 0.5. We set ri = T2 — 1/2 and g = ga 
(symmetric cloning). The dashed line is the value 2/3. The 
fidelity does not depend on the displacement amplitude a. 

This fidelity is plotted in Fig. 13 as a function of the 
squeezing parameter for different values of rj. We clearly 
see that for coherent states (corresponding to ^=0), the 
fideHty is 2/3 while decreasing with the degree of squeez- 
ing, eventually reaching zero for highly squeezed input 
states. 

In order to calculate the average fidelity, we assume 
that the squeezed state is drawn from an ensemble 

of states with a priori probability p{a,S,) = p{a)p{^). 
Above we mentioned that the cloning action with unity 
gain is independent on the distribution p{a), which can 
then be left undefined. The distribution of the squeezing 
factor is however quite important: it is clear that for 
completely unknown input squeezing (corresponding to 
a flat distribution) the average fidelity goes to zero. We 
therefore must restrict the set of input squeezed states 
to, say, a Gaussian distribution given by 

ft(0 = Aexp(-^j . (47) 

As evident from this expression we assume the distri- 
bution to be centered at ^ = which corresponds to a 
coherent state. This means that the coherent state is the 
most likely member in the set of input states, and we 
therefore conjecture that our machine is optimal in the 
Gaussian scenario. If however the distribution is centered 
at a known squeezing amplitude, say ^ = Coj then we 
believe that the optimal machine is the one mentioned 
above where the input states are unsqueezed [<S'~^(^o)] 
before the cloning machine and squeezed [S'(fo)] again 
after the cloning action. 

Using the polar coordinates, = pdpd(f), ^ = pe*"^, 
and Fjf{£^) — -FijdCI), the average fidelity now reads 

F.,^ [ d'^piOF.iO (48) 
Jc 

= 2 I dp^e^vi-^] F^{p) (49) 
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FIG. 4: Plot of the average fidelity of a set of squeezed 
states as a function of (see text for details) and different 
values of the efficiency rj: form top to bottom 77 = 1.0, 0.75, 
and 0.5. The dashed line corresponds to 2/3, i.e., the optimal 
cloning fidelity of coherent states. We put t\ = T2 = 1/2 and 
g = fls- 



This function is depicted in Fig. ^ as a function of for 
different values of ?/. If the standard deviation Ug ~ 0, 
the distribution in (|47|l is a delta function and the input 
alphabet contains only coherent states. In this case it 
reduces to the case discussed in the previous section and 
the expected fidelity is 2/3 (for ideal detection efficiency) 
as seen in the figure. We also see that the fidelity de- 
grades as the width of the distribution of the squeezing 
parameter increases, and eventually reaches zero when 
the a priori information is poor. At this point we should 
note that if one allows for non-Gaussian output clones 
the fidelity can be improved. E.g. it is known that the 
optimal doner of coherent states and the optimal uni- 
versal doner employ non-Gaussian operations and they 
yield fidelities of 68.3% ,2J] and 50% respectively. 



D. Thermal states 

Another interesting class of Gaussian states is the set of 
displaced thermal states gth.a = D[a)vt)xD'^{c^)^ which 
arise, for example, from the propagation of coherent 
states in a noisy environment [25|- The thermal state 
fth is given by 

■m— ^ ^ 

where N is the average number of thermal photons. Its 
covariance matrix is given by (Tin = {N + Since 
vth and, in turn, D{a) vth -D^(a) are not pure states, the 
cloning fidelity (TV) should be calculated using the full 
expression of Eq. H3t)|l . and the result is plotted in Fig. [S] 
as a function of TV and different values of 77. For the unity 
gain doner and assuming the detection efficiency to be 
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FIG. 5: Plot of the input-output fidelity F,,{N) of the dis- 
placed thermal state Qth,a as a function of the average number 
of thermal photons A'' and different values of the efficiency ry: 
form top to bottom rj — 1.0, 0.75, and 0.5. The dashed line 
corresponds to 2/3, i.e., the optimal cloning of coherent states. 
We put Ti — T2 = 1/2 and g — Qs. 
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FIG. 7: Plot of the average fidelity Fn of the set of thermal 
states distributed according to a "half- Gaussian" distribution 
l)54|l as a function of /ijv and different values of the efficiency 
T): form top to bottom r] = 1.0, 0.75, and 0.5. The dashed 
line corresponds to 2/3. We put t\ = T2 = 1/2 and g — gs- 
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FIG. 6: Plot of the average fidelity Fn of the set of thermal 
states distributed according to the top-hat distribution l|53|l 
as a function of the threshold value TV and different values of 
the efficiency rj: form top to bottom rj = 1.0, 0.75, and 0.5. 
The dashed line corresponds to 2/3. We put ti = T2 — 1/2 
and g = gs- 



where 

I otherwise 
for a top-hat distribution, and 

2 ( ^ 

V27r^^ L 2^A, J 

for a (re-normalized) "half-Gaussian" distribution. In 
Figs. Eland [3 we show the corresponding average fideU- 
ties, as functions of A/" and fiN, respectively, for different 
values of rj. For the top-hat distribution the average fi- 
delity monotonically increases as the threshold value A/" 
increases, whereas for the half-Gaussian one the average 
fidelity shows a maximum value depending on the value 
of r], as far as ?7 > 0.7. 



ideal (77 = 1), we derive the expression 

F,=i(iV)=Q + iV(3 + 2Ar) 

- ^/N{2N + l)(2iV2 + 5N + 3)^ . (51) 

We see that the fidelity increases with the average num- 
ber of thermal photons, that is, using the fidelity as a 
measure, the quality of the cloning action increases with 
the mixedness of the input states. 

Let us now consider a different ensemble of displaced 
thermal states, with random displacement and average 
number of thermal photons N distributed around zero 
either as a bounded flat, top- hat, distribution or as a 
"half-Gaussian" distribution. The average fidelity is 

F^^ / dNp{N)F^{N), (52) 
Jo 



IV. CONCLUSIONS 

We have analyzed in details a recently demonstrated 
scheme for linear cloning of Gaussian states |1J|. Using 
a suitable phase-space analysis the input-output fidelity 
has been evaluated for a generic (pure or mixed) Gaus- 
sian state taking into account the effect of non-unit quan- 
tum efficiency of homodyne detection and fluctuations in 
the beam splitters transmittivity. Our results indicate 
that the linear cloning machine suggested in is ro- 
bust against fluctuations of transmissivity and non-unit 
quantum efhciency 

We have explicitly evaluated the cloning fidelity for 
specific classes of non coherent displaced states. We 
found that a fixed (unknown) squeezing of the input 
states degrades the fidelity with respect to the coherent 
level, as one may expect for cloning of highly nonclassical 
states, while, on the contrary, cloning of displaced ther- 
mal states may be achieved with larger fidelity. Using 
the above results we have evaluated the average cloning 
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fidelity for classes of Gaussian states with fluctuating co- 
variance matrix, as for example displaced squeezed or 
displaced thermal states with the degree of squeezing or 
the number of thermal photons randomly distributed ac- 
cording to a Gaussian or a uniform distribution. Results 
indicate that the average fidelity monotonically decreases 
as the squeezing dispersion increases, whereas the be- 
haviour with respect to dispersion of thermal photons is 
not monotone. 
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